电导率各向异性的海洋电磁三维有限单元法正演   总被引：10，自引：8，他引：2       下载免费PDF全文

蔡红柱  熊彬  Michael Zhdanov 《地球物理学报》2015,58(8):2839-2850

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.  相似文献   

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany          下载免费PDF全文

Jana H. Börner  Feiyan Wang  Julia Weißflog  Matthias Bär  Ines Görz  Klaus Spitzer 《Geophysical Prospecting》2015,63(6):1430-1449

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

Dynamic FE simulation of four‐story steel frame modeled by solid elements and its validation using results of full‐scale shake‐table test          下载免费PDF全文

Tomoshi Miyamura  Takuzo Yamashita  Hiroshi Akiba  Makoto Ohsaki 《地震工程与结构动力学》2015,44(9):1449-1469

Dynamic finite element analyses of a four‐story steel building frame modeled as a fine mesh of solid elements are performed using E‐Simulator, which is a parallel finite element analysis software package for precisely simulating collapse behaviors of civil and building structures. E‐Simulator is under development at the National Research Institute for Earth Science and Disaster Prevention (NIED), Japan. A full‐scale shake‐table test for a four‐story frame was conducted using E‐Defense at NIED, which is the largest shaking table in the world. A mesh of the entire structure of a four‐story frame with approximately 19 million degrees of freedom is constructed using solid elements. The density of the mesh is determined by referring to the results of elastic–plastic buckling analyses of a column of the frame using meshes of different densities. Therefore, the analysis model of the frame is well verified. Seismic response analyses under 60, 100, and 115% excitations of the JR Takatori record of the 1995 Hyogoken‐Nanbu earthquake are performed. Note that the simulation does not reproduce the collapse under the 100% excitation of the Takatori record in the E‐Defense test. Therefore, simulations for the 115% case are also performed. The results obtained by E‐Simulator are compared with those obtained by the E‐Defense full‐scale test in order to validate the results obtained by E‐Simulator. The shear forces and interstory drift angles of the first story obtained by the simulation and the test are in good agreement. Both the response of the entire frame and the local deformation as a result of elastic–plastic buckling are simulated simultaneously using E‐Simulator. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Finite‐element analysis of laminated rubber bearing of building frame under seismic excitation          下载免费PDF全文

Makoto Ohsaki  Tomoshi Miyamura  Masayuki Kohiyama  Takuzo Yamashita  Masashi Yamamoto  Naohiro Nakamura 《地震工程与结构动力学》2015,44(11):1881-1898

Finite element analysis is carried out for a building frame supported by laminated rubber bearings to simultaneously investigate global displacement and local stress responses under seismic excitation. The frame members and the rubber bearings are discretized into hexahedral solid elements with more than 3 million degrees of freedom. The material property of rubber is represented by the Ogden model, and the frame is assumed to remain in elastic range. It is shown that the time histories of non‐uniform stress distribution and rocking behavior of the rubber bearings under a frame subjected to seismic excitation can be successfully evaluated, and detailed responses of base and frame can be evaluated through large‐scale finite element analysis. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

应力扰动对地震周期和地震矩影响的数值模拟          下载免费PDF全文

翁辉辉  黄金水 《地震学报》2015,37(1):65-79

利用二维有限元数值模型,结合断层滑移弱化摩擦准则对断层滑动规律以及应力扰动对其影响进行了研究．数值计算结果表明,在均匀应力分布情况下, 平面断层滑动显示出典型的特征地震规律,断层面上的应力扰动对断层滑动规律产生影响,压应力增加明显延迟地震的发生时间,并增加地震释放的能量．应力扰动发生在地震破裂临界区时的影响比在震前滑移区时的影响显著．当发生在地震滑移区时,若应力扰动足够大,则压应力增大会造成地震发生时部分动力断层被暂时锁住,使得地震释放的能量变小,但可增加后续地震的能量; 而压应力减小则可导致地震规律产生更加复杂的变化,会即时触发地震．如果应力扰动发生在一个地震周期的早期,则触发的地震较小,但可导致随后的地震提前发生; 如果应力扰动发生在一个地震周期的后期,则会触发大地震．当应力扰动位于震前滑移区或破裂临界区时,小的扰动也可能产生类似的效果．应力扰动产生越晚,这种影响也越明显．应力扰动发生在破裂临界区的影响最明显．应力扰动的影响一般主要集中在应力发生扰动后的1—2个地震周期内．后续地震基本恢复无应力扰动时的特征地震规律．   相似文献   

three-dimensional sensitivity coefficients of apparent resistivity and preliminary application     

XIE Tao  LU Jun 《地震地质》2015,37(4):1125-1135

We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.  相似文献   

金融空间联系及K-means聚类中心等级识别研究——以长三角为例     

杨志民  化祥雨  叶娅芬  邵元海 《地理科学》2015,35(2):144-150

以2001年、2006年、2011年长三角城市金融机构人民币存款、贷款额数为样本,构建金融空间联系模型,定量分析长三角城市金融空间联系分异特征。在此基础上构建K-means金融中心等级识别模型,识别长三角城市金融中心等级。研究表明：① 2001-2011年长三角城市金融“质量”空间趋势较为稳定,总体呈现东部高于西部,中部高于南、北部的倒U形分布。② 金融空间联系最大引力线联结格局较为稳定。③ 金融空间联系网络结构格局变化显著,主要从简单的“折线型”空间结构逐渐发展成简单的、复杂的“网络型”空间结构。④ 长三角金融城市中心等级空间分布格局稳定,以上海市金融中心最为突出。  相似文献   

基于要素流的长三角城市群空间布局     

顾雯娟  欧向军  叶磊  杨宝宝 《热带地理》2015,35(6):833-841

以长三角城市群30个地级市市区为研究单元,从要素流动的动态视角,建立要素流的数理模型,借助于ArcGIS软件,对长三角城市群内城市的经济流、交通流和信息流的强度和方向进行定量分析和评价,运用赋值法对3种要素流总量进行评分,并进行空间层级划分。得出核心区（包括上海、江苏沿江八市、环杭州湾地区等地市）和辐射影响区（位于长三角两翼的地市）。基于城市间的空间联系特征和城市群的空间层级结构,对长三角城市群的空间发展布局提出进一步引导和规划,即上海为核心,苏州、南京、杭州、合肥和徐州为区域性中心城市,构成以沪宁合、沪杭、沪徐和沪温为4个发展轴的反“K”字型的空间形态,依据城市区位和要素联系特征,将城市群进一步划分为6个城市组团,形成以中心城市为增长节点、以发展轴为联结线、以城市组团为发展域的“点―线―面”有序发展的空间互动格局。  相似文献   

铅垂断层向错引起的地表重力变化     

段虎荣  杨凡 《大地测量与地球动力学》2015,35(4):627-631

利用直立长方体组合模型以及坐标旋转理论,对铅垂断层向错引起的重力变化进行数值模拟,分别计算深度为1、5、10 km的铅垂断层在单一方向上、两个方向上向错引起的地面重力变化。结果表明,当断层在单一方向产生向错时,以左端点为基点的W1、W2向错引起的地面重力变化场均呈现以断层为轴线的对称分布|以中点为基点的W3向错引起的地面重力变化场构成四象限反对称分布。当断层在两个方向上同时产生向错时,以左端点为基点的W1W2向错引起的地面重力变化场呈现以断层为轴线的对称分布,整体显示W2向错起主导作用|以中点为基点的W2W3向错引起的地面重力变化场呈现非对称分布,整体也显示出W2向错起主导作用,断层附近的扭错现象为W3向错的作用。随着断层深度的增加,断层向错引起的地面重力变化值递减,但影响范围逐渐扩大。  相似文献   

北斗三频无几何、消电离层组合周跳探测方法研究     

王华润  柴洪洲  谢恺 《大地测量与地球动力学》2015,35(3):406-411

根据多频组合理论,构造基于北斗三频数据的无几何、消电离层组合,结合双频MW组合法和电离层残差法,对北斗实测三频数据进行周跳探测。实验表明,上述3种组合的结合使用能探测出所有大小的周跳,并具有较好的修复效果。  相似文献